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ABSTRACT 

Context. Photometric variability is a distinctive feature of young stellar objects; exploring variability signatures at different wave¬ 
lengths provides insight into the physical processes at work in these sources. 

Aims. We explore the variability signatures at ultraviolet (UV) and optical wavelengths for several hundred accreting and non¬ 
accreting members of the star-forming region NGC 2264 (~3 Myr). 

Methods. We performed simultaneous monitoring of u- and r-band variability for the cluster population with CFHT/MegaCam. The 
survey extended over two full weeks, with several flux measurements per observing night. A sample of about 750 young stars is 
probed in our study, homogeneously calibrated and reduced, with internally consistently derived stellar parameters. Objects span the 
mass range 0.1-2 Mq; about 40% of them show evidence for active accretion based on various diagnostics (H^, UV, and IR excesses). 
Results. Statistically distinct variability properties are observed for accreting and non-accreting cluster members. The accretors 
exhibit a significantly higher level of variability than the non-accretors, in the optical and especially in the UV. The amount of u- 
band variability is found to correlate statistically with the median amount of UV excess in disk-bearing objects, which suggests that 
mass accretion and star-disk interaction are the main sources of variability in the u band. Spot models are applied to account for the 
amplitudes of variability of accreting and non-accreting members, which yields different results for each group. Cool magnetic spots, 
several hundred degrees colder than the stellar photosphere and covering from 5 to 30% of the stellar surface, appear to be the leading 
factor of variability for the non-accreting stars. In contrast, accretion spots with a temperature a few thousand degrees higher than the 
photospheric temperature and that extend over a few percent of the stellar surface best reproduce the variability of accreting objects. 
The color behavior is also found to be different between accreting and non-accreting stars. While objects commonly become redder 
when fainter, typical amplitudes of variability for accreting members rapidly increase from the r to the u band, which indicates a much 
stronger contrast at short wavelengths; a lower color dependence in the photometric amplitudes is instead measured for diskless stars. 
Einally, we compare the u-band variability monitored here on two-week timescales with that measured on both shorter (hours) and 
longer (years) timescales. We find that variability on timescales of hours is typically ~10% of the peak-to-peak variability on day 
timescales, while longer term variability on a timescale of years is consistent with amplitudes measured over weeks. 

Conclusions. We conclude that for both accreting and non-accreting stars, the mid-term rotational modulation by hot and cold spots 
is the leading timescale for a variability of up to several years. In turn, this suggests that the accretion process is essentially stable over 
years, although it exhibits low-level shorter term variations in single-accretion events. 

Key words. Accretion, accretion disks - stars: low-mass - stars: pre-main sequence - stars: variables: T Tauri - open clusters and 
associations: individual: NGC 2264 - ultraviolet: stars 


1. Introduction 

The notion of photometr ic variabili ty is at the very base of the 
T Tauri class definition (lJovlll94^ . Since the first pioneering 


* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-Erance-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRG) of Canada, the Institut National des Sciences de I’Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. 


work (see lHerbig|[i962h . variability appears to be a ubiquitous 
property of these young stellar objects (YSOs), as universal a 
feature as bewildering in its intrinsic case-to-case diversity. 

The variable nature of young stars manifests on a 
wide range of waveleng ths and on different timescales (see 
iMenard & Bertou3ll999l for a review). Intense X-ray emission 
and variability are exhibited by young cluster members both 
on short (~ hours) an d mid-term (~davs ) timescales, ensuing 
from magnetic flaring (iFavata et al.ll2005]; also seen in the UV) 
and/or modulation effects linked with rotation (iFlaccomio et al.l 
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1201(1 12012h . respectively. YSOs are found to be variable 
in th e optical on timesca les ranging from fractions of days 
(e.g., RucinskieL^|200 ^ to weeks, months, and years (e.g., 
iBouvier et ahllQQsC 1995UGrankin et al.l2007ll20()^ : these vari¬ 
ations reflect the geometry of the systems as well as the physics 
and intrinsic timescales pertaining to their dynamics. Variability 


on an assortment of timescales is also dete cted in the near- (e.g ., 
IWolk et ^l2013h and mid-infrared (e.g.. iGiinther et aP 120141) . 


where a major contribution to the measured flux arises from ther¬ 
mal emission from the inner disk, when present. 

An overall similarity in variability pattern appears to be 
shared by weak-lined T Tauri stars (WTTS), more evolved YSOs 
with no evidence of circumstellar disk; their light curves typi¬ 
cally show a regular, periodic profile, with amplitudes ranging 
from < 0.1 to ~ 0.5 mag in the optical, and are stable over tens 
to hundreds of rotational cycles. A more complex and diverse 
picture is observed for classical T Tauri stars (CTTS), whose dy¬ 
namics are dominated by the interaction with an active accretion 
disk. Optical light curves typically show large (up to a few mag) 
variability amplitudes and often irregular profiles, with semiperi- 
odic components coexisting with rapid, stochastic flux variations 
and fading events. As a common rule, TTS variability amplitudes 
are found to decrease toward longer wavelengths, albeit with dif¬ 
ferent rates for different types of sources. 

A first thorough exploration of optical YSO variability and 
of its physical in terpretation was provided by the study of 
iHerbst et'^(ll994ll . Based on t/BTf?/monitoring of 80 sources, 
covering tens of epochs over several years, the authors catego¬ 
rized three main types of mid-term TTS variability. Type I, re¬ 
sponsible for the well-behaved, periodic light curves observed 
for WTTS, consists of a simple modulation effect produced by 
cool magnetic spots at the stellar surface during stellar rotation. 
Except for episodes of magnetic flaring, this is sufficient to ex¬ 
plain the variability of disk-free YSOs. This also represents an 
underlying component in the variability of their similarly mag¬ 
netically active, disk-bearing counterparts, where the physics of 
accretion and star-disk interaction plays a prominent role. Type 
II variability, specific to CTTS, originates in a changing mixture 
of cool magnetic spots and hot accretion spots at the surface of 
the star. Hot spots are the signatures of the impact of the accre¬ 
tion column from the inner disk onto the central object; the en¬ 
ergetic emission produced in the shock produces a flux excess at 
short wavelengths (UV, soft X-rays), and hence these sources are 
bluer than photospheric colors would suggest. Periodic behavior 
can be observed in some cases, but irregular variation patterns 
are more common, which suggests shorter lifetimes for accretion 
spots than for magnetic spots. Type III variability is of a differ¬ 
ent nature, pertaining again only to disk-bearing sources. Objects 
in this category show light curves with typical brightness levels 
close to the maximum luminosity state, interspersed by dips of 
up to a few magnitudes, with little color dependence in the am¬ 
plitude. Variable obscuration from circumstellar dust appears to 
be the dominant source of variability. A quasi-periodicity in the 
fading events may be observed (e.g.. iBouvier et aHl2007h : this 
is indicative of a concentration of material in the inner disk that 
rotates with the system and periodically occults part of the pho¬ 
tosphere. 

From this discussion, it is clear that large-scale, synoptic in¬ 
vestigations of YSO variability signatures over the full wave¬ 
length range are critical for understanding the different nature of 
these systems. In addition, a time sampling on timescales short 
enough to trace the short-term variations and extending over 
timescales relevant to various processes is needed to achieve a 
detailed physical description of the processes at play at the stel¬ 


lar surface and in the star-disk interface. Remarkable progress in 
this respect was marked by the advent of space-based telescopes, 
with their exquisite cadence and photometric accuracy. The po¬ 
tential of s uch campaigns was illustrated by the Spitzer Space 
Telescope (IWerner et alJl2004l) InfraRe d Array Camera (IIT4C; 


Fazio et ^ 2004) YS(3 VAR project (iMorales-Calderon et al.l 


201 iblRebull et al.ll2014ll and by a first investigation of the star¬ 
forming region NGC 2264 (~3 Myr) with the COnvection, RO- 
tation and planetary Transits ( CoRoT) optical space telescope 
(lAlencar et al.ll2010l) . In December2011, the improved capabil¬ 
ity of space-based surveys in time-domain exploration and the 
advantages of a multiwavelength approach to YSOs variability 
were combined in the unprecedented “Coordinated S ynoptic In¬ 
vestig ation of NGC 2264” (CSI 2264) campaign (ICodv et al.l 
l2014h . This consisted of a unique cooperative observing pro¬ 
gram aimed at mapping the photometric variability of the whole 
cluster (~ 1500 known members) simultaneously in infrared, op¬ 
tical, UV, and X-ray wavelengths on timescales varying from 
hours to over a month. The primary components of the opti¬ 
cal and infrared observations were performed with CoRoT and 
Spitzer/IRAC for a duration of 40 and 30 co nsecutive days, re - 
spectiyely. This dataset, described in detail in ICody et al.l (l2014ll . 
enabled a more accurate classification of the diyerse nature of 
disk-bearing objects based on light-curye morphology, where de¬ 
fined classes included quasi-periodic, dippers, bursters, stochas¬ 
tic behayior, strictly periodic, and long-timescale yariations. 
Complementary information gathered at similar and at different 
wayelengths with 13 additional space-borne and ground-based 
telescopes during the campaign allowed us to inyestigate the 
yariable nature of these objects in more detail. 


The UV wayelength domain offers a particularly interesting 
perspectiye on the star-disk interaction dynamics, as it proyides 
one of the most direct diagnostics of disk accretion onto the 
central object. The association between UV excess displayed 
by TTS and infall of material onto the star from a surround- 
ing reg i on of moderate yertical extent was first suggested by 
IWalken (Il972h . about twenty years before magnetospheric ac- 
cretion models for dis k accretion in CTTS were deyeloped (e.g., 
iHartmann et al.l[T994ll . In the current picture, accretion columns 
impact the stellar surface at near free-fall yelocities, hence result¬ 
ing in shock regions (hot spots) close to the magnetic poles. The 
energetic emission produced in these areas results in the distinc- 
tiye UV flux excess detected for accreting stars compared to their 
non-accreting counterparts. Remarkably, yery little systematic 
inyestigation of YSO yariability in the UV exists in the literature 
to date, particularly studies encompassing large statistical popu¬ 
lations of different nature with time coyerage and sampling suit¬ 
able to qualify different component s. An interesting effort in thi s 
respect is proyided by the study of iFallscheer & HerbstI (l2006l) . 
who monitored the [/-band yariability of about a hundred young 
stars in the NGC 2264 held to inyestigate disk-locking, using 
the UV excess as a proxy for accretion. Within the context of 
the CSI 2264 project, we performed a dedicated UV monitor¬ 
ing campaign of seyeral hundred young members of NGC 2264 
in February 2012, simultaneously with optical monitoring, at 
the Canada-France-Hawaii Telescope (CFHT). For the first time, 
this suryey proyides direct measurements of mass accretion rates 
(Macc) and of their yariability across a large statistical sample of 
coeyal pre-main sequence (PMS) stars, all similarly character¬ 
ized. The extent and homogeneity of the dataset allows us to i) 
inyestigate the dependence of the mass accretion rate on stellar 
parameters such as mass and eyolutionary age, ii) explore the 
nature of the wide spread in Macc yalues detected in any mass 
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bin, iii) evaluate the typical amount of variability on Macc on 
timescales relevant to stellar rotation, iv) assess the contribution 
of geometric modulation effects to the detected am ount of M„rr 
variab ility. Results of this analysis were reported in I Venuti et al.l 
(1201 4 hereafter Paper I). 

In a complementary effort, here we focus on exploring and 
characterizing the variability signatures at short wavelengths that 
pertain to different types of YSO variables. The paper is orga¬ 
nized as follows. Section 12 provides a description of the mon¬ 
itoring survey and of the subsequent light-curve extraction and 
processing. Section irTI presents a comparison of the optical and 
UV variability level exhibited by accreting and non-accreting 
sources, which were probed using different variability indices. 
Variability and accretion indicators are combined in Sect. l3.2l 
to evaluate the impact of the accretion process on the substan¬ 
tial variability characteristic of CTTS. In Sect. 13.31 we explore 
the color signatures associated with different physical processes 
at the origin of YSO variability, that is, how steeply the vari¬ 
ability amplitudes change from the optical to the UV. We show 
detailed color-magnitude-correlated variations for representative 
sources of different variability types and highlight how these di¬ 
agrams reflect the underlying physics. In Sect. 13.41 we adopt a 
spot model description of the observed optical -i- UV variability 
amplitudes to derive a statistical depiction of the distinctive sur¬ 
face features of accreting vs. non-accreting objects. Our results 
are summarized and discussed in Sect.H] while Sect. |5] presents 
our conclusions. 

2. Observations and data analysis 

We extensively monitored the u- and r-band variability of 
NGC 2264 members from February 14 to 28, 2012, at the CFHT 
using the wide-held camera MegaCam (0.96° x 0.94° FOV). A 
full description of the monitoring survey and of the subsequent 
data processing and light-curve extraction has been provided in 
Paper I. Here, we summarize the main features of the observing 
campaign; we refer to Paper I for further details. 

Targets were observed during 11 nights distributed along the 
14-day-long run; among these, 7 nights were photometric. A sin¬ 
gle u + r observing block, consisting of hve dithered exposures in 
u and hve dithered exposures in r, obtained consecutively, was 
performed repeatedly on each observing night during the run; 
each exposure was then individually processed to reconstruct the 
light variations of all sources on different timescales (hours to 
days). 

The large population of held stars detected in the telescope 
aperture (~8000 objects) was used to test the accuracy of the 
photometry of individual exposures by comparison with the mas¬ 
ter frame; photometric sequences obtained under poorer observ¬ 
ing conditions (hence resulting in a large magnitude offset from 
the master frame) were rejected upstream. A log of the Mega- 
Cam monitoring survey of NGC 2264 is provided in Table [T] 

A 5(T-clipping procedure was adopted to identify and re¬ 
ject occasional very discrepant points in the light curves; this 
occurred in less than 25% of the cases in the M-band and for just 
~5% of the r-band light curves. 

The catalog of NGC 2264 members, assembled as described 
in Paper I, was used as the reference for the identihcation of 
members and held stars in the monitoring survey. Instrumen¬ 
tal photometry was statistically calibrated to the Sloan Digi¬ 
tal Sky Survey (SDSS) photometric system, using a popula¬ 
tion of -3000 stars in the CFHT FOV for which SDSS pho- 
tomet ry is publicly available from the seventh SDSS Data Re¬ 
lease (lAbazaiian et al.ll2009h . Our census of NGC 2264 mem- 


Table 1. Log of CFHT/MegaCam monitoring survey of NGC 2264 
(Feb. 14-28, 2012). 


Date 

Flagpho/ 

n° obs. seq/ 

phot 




12.1 

Feb. 15 

1 

3 

72.8 




72.9 




73.7 

Feb. 16 

1 

3 

73.8 




73.9 

Feb. 17 

0 

4 


Feb. 18 

0 

7 


Feb. 19 

0 

3 





78.7 

Feb. 21 

1 

3 

78.8 




78.9 




81.7 

Feb. 24 

1 

3 

81.8 




81.9 

Feb. 26 

2 

4 

83.7 

83.8 

Feb. 27 

2 

4 

84.8 

84.9 

Feb. 28 

2 

3 

85.7 

85.8 

Feb. 29 

0 

1 



Notes. Flag on observing conditions: 1 = night was photometric; 2 = 
parts of the night were photometric; 0 = night was non-photometric. 

Number of observing sequences performed during the night (i.e., 
number of times the same observing block was repeated during the 
night). Julian Date (JD — 2 455 900, truncated at the Hrst decimal 
place) of observing sequences obtained in photometric conditions dur¬ 
ing the night. 


bers is complete down to « ~ 21.5 and r ~ 18.5 mag, with 
the brightest objects reaching apparent magnitudes of m ~ 15 
and r ~ 13.5. A population of ~750 members has been probed, 
spanning 0.1-2 Mq; of these, about 40% are actively accreting. 
Criteria adopted for member classification are detailed in Paper I 
and mainly include H„ emission, variability, and IR excess from 
both our and earlier surveys; for this part of the analysis, UV 
excess information from our data was used as a complement to 
refine the census of NGC 2264 members with the identification 
of previously unknown accreting sources in the cluster. Here¬ 
after the objects with accretion disk signatures are referred to as 
CTTS, while those lacking any of them are called WTTS. A full 
list of stellar parameters for the targeted objects is provided in 
Table 3 of Paper I; CTTS/WTTS classification and spectral type 
estimates are reported in Table |2of this paper. 

Two main timescales contribute to the monitored variabil¬ 
ity: i) short-term (< hours), dominated by short-lived events (like 
flares or bursts), and ii) mid-term (—rotation period), dominated 
by geometric modulation effects. The latter is the component of 
interest for the present study. To smooth out short-term variabil¬ 
ity and the photometric noise components from single observ¬ 
ing sequences, we computed the 10'^' and 95^^ percentile lev¬ 
els, in magnitudes, over the whole light curve and disregarded 
all points outside this range. These selective levels were set em¬ 
pirically upon examination of the full statistical sample of light 
curves; a more severe threshold was adopted at the brighter end 
to properly omit flares. 
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Fig. 1. Two examples of CFHT r-band light curves at the various steps 
of processing (detailed in Sect. [He gray points indicate observations ob¬ 
tained in non-photometric nights or sequences, which were rejected for 
the present analysis; light-curve points excluded after 10'*-95'* per¬ 
centile selection are marked in red; the final set of data points retained 
is shown in blue. 

The various steps of light-curve processing are illustrated in 
Fig.dJfor two cases. Similar diagrams for individual objects were 
visually inspected at the end of the routine to verify accurate 
point selection. In the few tens of cases for which incorrect au¬ 
tomatic selection resulted in biased amplitudes (e.g., when light- 
curve minima or maxima occurred close to non-photometric 
nights), this was corrected by hand. 

3. Results 

3.1. u-band variability of young stars in NGC 2264 

Measuring the light-curve dispersion around the average pho¬ 
tometric level provides a first snapshot of variability properties 
across the sample. Non-variable stars are expected to have weak 
stochastic fluctuations around their mean brightness due to the 
photometric noise that affects every measurement. On the other 
hand, light curves of intrinsically variable stars are expected to 
exhibit a significantly larger root-mean-square (rms) variation 
than that expected for photometric noise. 

These populations of non-variable and variable stars can be 
observed in Fig. [2] for main-sequence dwarfs and young mem¬ 
bers in the NGC 2264 held. The light curve rms measured for 
field stars as a function of magnitude defines both in r and in 
u the photometric noise level at a given brightness. In each (m, 
rms„,) diagram, the young stars exhibit a distinct level of intrin¬ 
sic variability over the photometric noise sequence; furthermore, 
accreting and non-accreting members tend to occupy statisti¬ 
cally distinct but overlapping regions, with a median rms value 
2.5 times higher for the first group than for the second, and in¬ 
creasing in both cases from the r to the u band. A two-sample 
Kolmogorov-Smirnov test, aimed at comparing the distribution 
in rms of accreting members to that of non-accreting objects both 


Fig. 2. r-band (left panel) and j(-band (right panel) light-curve disper¬ 
sion as a function of magnitude for different groups of monitored ob¬ 
jects: field stars (black crosses), accreting NGC 2264 members (CTTS; 
magenta circles), and non-accreting members (WTTS; blue squares). 

in u and in r, attests to the significantly different nature of the 
two populations, with a null-hypothesis probability of 10“^^ in 
r and 10“^^ in u. For fainter objects, the photometric noise af¬ 
fecting the light-curve measurements becomes more important, 
especially in the u band, which results in a certain degree of over¬ 
lapping between the loci of variable and non-variable stars. 

A more accurate variability tracer is provided bv IStetsonI ’s 
(1 19961) J index. This probes the level of correlated variability be¬ 
tween different wavelengths monitored at the same time, which 
significantly reduces the spurious contribution from stochastic 
noise to the detected variability level. Points in the u-band light 
curve are paired with the closest r-band epoch (at an average 
separation of 10 to 15 minutes). The J index is then defined as 

, SLisgn(P/) 

L,=i 

where i is the time index across the series, n is the total number 
of simultaneous u, r pairs (typically 98 over t wo weeks), “sgn” i s 
the sign function, ic, is a weight, defined as in iFruth et al.l(l2012h . 
to account for the actual time distance between paired u and r 
measurements, and P,- is the product of normalized residuals of 
paired observations: 


” 1 

f Hi M (J- \ 

/ r,— < r >o.\ 

n — 1 

V o-,(m) j 

V o-iir) J 


(m, and r, being the magnitudes at time i, < u >o- and < r >o- 
the mean magnitudes and cr the photometric uncertainty). 

A spurious, noise-driven “variability” component is expected 
to statistically average to zero throughout the whole time series. 
Conversely, a definite correlation is expected to exist between u- 
and r-band behavior for a truly variable star; this would result in 
values of J significantly different from zero. 

Figure|2 shows the J distribution for cluster members and 
field stars in our survey. Field stars are predominantly located on 
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Fig. 3. IStetsonl ’s (Il996h J index of correlated UV/optical variability as 
a function of n-band magnitude for different groups of monitored ob¬ 
jects: field stars (black crosses), accreting NGC 2264 members (CTTS; 
magenta circles), and non-accreting members (WTTS; blue squares). 


Fig. 4. UV excess measurements (see Paper I) are compared to the 
level of M-band rms variability detected for CTTS (magenta circles) and 
WTTS (blue squares) in NGC 2264. Objects displayed are all selected 
as significantly variable using the J-index indicator (see text). The UV 
excess increases from right to left on the diagram. 


a narrow, horizontal sequence centered on 7=0, independent of 
the object’s brightness. A low percentage of field stars (~3%) 
shows a significant amount of variability, attributable either to 
occasional spurious points or to an intrinsic variable nature (e.g., 
binaries, pulsators, active stars). Of the NGC 2264 members, 
about 63% show clear variability above the 7 sequence traced by 
field stars at the 3 cr level (75% at 1 cr level); this percentage rises 
to 72% (81%) when considering the accreting members alone. 
As for Fig. 121 CTTS display a significantly higher level of vari¬ 
ability than WTTS, with the typical 7 value three to four times 
higher than that of non-accreting members. 

A variability flag, median u- and r-band photometry, light- 
curve amplitudes, light curve rms, and 7 index values for indi¬ 
vidual members are reported in Tabled 

3.2. UV variability and accretion 

As discussed in Sect. 13. II CTTS are statistically found to exhibit 
significantly stronger variability than WTTS; this is especially 
observed at shorter wavelengths (m band). This suggests that ac¬ 
cretion mechanisms, whose signatures can be most directly de¬ 
tected in the UV, are a primary cause of the enhanced variability 
of these young stars. We would then expect to observe a direct 
link between accretion diagnostics and variability indicators. 

In Paper I, we reported on a full characterization of accretion 
properties for NGC 2264 members from the UV excess diagnos¬ 
tics; we showed that the UV color excess displayed by accreting 
stars relative to WTTS provides a direct proxy for the accretion 
luminosity. 

In Fig.|4] we compare UV excess measurements to the level 
of M-band rms variability. The diagram is obtained as follows: 

(a) We selected from Fig.[3]all CTTS (217) and WTTS (210) 
whose J index was at least 3 cr higher than the typical 
location at the corresponding brightness on the field star 


sequence. For a given object, the measured light curve rms 
is the result of two independent contributions, intrinsic vari¬ 
ability and photometric noise; cr„ea., = ^rrltj + rrloise ■ 
recover rrobj, we estimated the lowest magnitude-dependent 
noise component expected to affect the measurement as a 
fit to the lower envelope of the field star locus in Fig.|2] and 
quadratically subtracted this contribution from the measured 
rms values for young stars. 

(b) UV excesses were measured as in Eq. 9 of Paper I for 
CTTS. The same procedure was applied here to compute 
a “UV excess” for WTTS, for comparison purposes. In 
WTTS, the UV excess definition has no direct meaning, 
but it provides a measure of the scatter of WTTS colors 
around the reference sequence. Enhanced chromospheric 
activity is expected to translate both into somewhat bluer 
colors and stronger variability. Recent studies conducted on 
some what older PMS po pulations than NGC 2264 members 
(e.g.. iKamai et al.ll20l4 on the Pleiades) have shown that 
more rapidly rotating objects might display bluer colors and 
stronger variability amplitudes in light curves. Hence, some 
correlation might be detected between color excess and 
amount of variability exhibited by WTTS. Eor the reasons 
discussed in Sect. 4 of Paper I, objects with r < 14.5 are not 
displayed in Pig. ID 

Sp earman’s non-parametric rank correlation test (iPress et al.l 
1 19921) was adopted to probe the statistical significance of any 
relationship between UV color excess and M-band rms variabil¬ 
ity, both in the accreting and non-accreting groups in Pig. ID A 
/7-value of 0.07 is returned from the test for the latter. More im¬ 
portantly, direct evidence of a correlation between the accretion 
process and photometric variability at short wavelengths is pro- 


Article number, page 5 ofll4l 


















A&A proofs: manuscript no. Venuti_UV_variability_NGC2264 


Table 2. Median photometry, variability amplitudes, light curve rms, J index, and color slopes for members discussed in this study. 


CSIMon-#‘ 

Var^ 

Class^ 

SpT^ 

medu 

medr 

ampu 

ampr 

nnsu 

nnsr 

J index 

Ar 

^ 5 

^ slope 

A(u-r) 

000007 

* 

c 

K7 

17.22 

14.49 

1.104 

0.207 

0.34 

0.05 

32.97 

0.174 

0.004 

000009 


w 

F5 

18.21 

14.99 

0.019 

0.010 

0.005 

0.003 

0.11 



000011 

* 

c 

K7 

16.76 

15.09 

1.008 

0.429 

0.34 

0.14 

58.26 

0.642 

0.008 

000014 


w 

K7:M0 

18.88 

14.25 

0.023 

0.007 

0.006 

0.003 

0.20 



000015 


w 

K7:M0 

19.39 

14.55 

0.028 

0.008 

0.008 

0.003 

0.60 



000017 

* 

c 

K5 

18.41 

14.94 

0.216 

0.130 

0.06 

0.04 

8.50 

1.52 

0.15 

000018 

* 

w 

K3:K4.5 

19.05 

14.99 

0.192 

0.073 

0.03 

0.02 

2.05 

0.70 

0.08 

000020 

* 

w 

K7 

20.02 

15.94 

0.153 

0.109 

0.04 

0.03 

1.65 

1.2 

0.3 

000022 


w 

M4 

21.80 

17.72 

0.238 

0.044 

0.06 

0.016 

0.20 



000023 


w 

M3 

22.30 

18.18 

0.475 

0.063 

0.11 

0.02 

0.07 




Notes. A full version of the Table is available in electronic f orm at the CDS. A portion is shown here for guidance regarding its form and content. 

Object identifiers adopted within the CSI2264 project fsee lCodv et al.l20l^ . Variability flag; objects that exhibit a distinct level of variability 
stronger than that for field stars, based on the J-in dex indicator , are m arked with an asterisk. “c” = CTTS; “w” = WTTS (classification from 
IVenuti et al.ll2014l) . Spectral type estimates from IVenuti et al.l ( l2014ll . ® rms uncertainty on Ar/A{u — r) color slope. 


vided by the distribution of accreting objects on this diagram. 
The bulk of the CTTS are located at a higher u-band rms level 
than the WTTS, as already observed in Fig.|2l and higher rms 
values correspond, on average, to higher (i.e., more negative) 
UV excesses. This correlation is estimated to be significant at 
the 6 cr level. 

3.3. UV variability and color signatures 

Figure|4]shows the global accretion-variability connection across 
the cluster. This statistical information, however, is not all that is 
needed to achieve a detailed picture of how variability relates 
to the physics of individual objects, as illustrated by the wide 
spread of rms (u) values at any given UV excess. 

More accurate information on the nature of YSO variabil¬ 
ity derives from investigating color variations. Indeed, while 
single-band amplitudes provide some indication on the contrast 
between stellar photosphere and the source of variability, color 
variations trace the wavelength dependence of this contrast and 
hence provide more direct clues to its physical origin. 

In Fig.|5] UV colors and variability properties for NGC 2264 
members in Fig.|4] are combined in a dynamic picture of the 
(m — r, r) diagram of the cluster. A significant number of CTTS 
display markedly bluer colors than the corresponding location on 
the cluster sequence traced by WTTS, as a result of the distinc¬ 
tive UV excess linked with ongoing disk accretion. In addition, 
as illustrated in Fig.H] high UV excesses are typically associ¬ 
ated with strong variability, whereas a significantly lower level 
of variability, both in magnitude and in color, is observed on the 
same timescales for non-accreting objects. 

The detailed picture of variability observed for individual ob¬ 
jects may vary on a case-by-case basis, as shown in Fig. 5. This 
is reflected in the range of slopes associated with variability bars 
among both WTTS and, especially, CTTS. Different color signa¬ 
tures are expected to correspond to different physical scenarios. 
When variability is driven by spot modulation, systems are ex¬ 
pected to be redder at fainter states, with variability amplitudes 
that decrease toward longer wavelengths; a steeper decrease with 
A is expected for hot spots than for cold spots. When variability 
is due to opaque material that crosses the line of sight to the star, 
little color dependence is expected in the main occultation event. 

This is well illustrated in the panels of Fig.0 which com¬ 
pares brightness and color variation trends for four members of 
NGC 2264 that are dominated by different mechanisms. Panel a) 



Fig. 5. Monitoring of color and magnitude variations on the (u — r, r) 
diagram for CTTS (magenta lines) and WTTS (blue lines) in NGC 2264 
on a timescale of weeks. Each bar represents a single object: the center 
of the bar is located at the color and r-band magnitude of the object at 
the median luminosity state, while the amplitudes along m —r and r mark 
the total amount of variahility detected during the CFHT monitoring 
(i.e., and over the 10"' -95"’ percentiles 

range, respectively). The effects of redden ing on the diagram are tra ced 
hy the black an'ow in the top right corner jFiorucci & Muna^l2003f) . 

shows an example of a stochastic accretor (IStauffer et al.ll201^ 
whose variability is dominated by short-lived, stochastic accre¬ 
tion bursts. The color light curve nearly retraces the features of 
the magnitude light curve, as the system becomes consistently 
bluer for increasing brightness. The same global property is ob¬ 
served for case b), which depicts a WTTS member whose vari¬ 
ability is produced by rotational modulation of cold magnetic 
spots. Cases b) and a) are distinguished by both i) the varia¬ 
tion amplitudes, especially in color, and ii) the slope of the av¬ 
erage r vs,, u — r trend (as can be deduced when comparing 
the direction traced by data points with the slope of the red- 
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a) Stochastic accretion bursts 


b) Cold spot modulation 
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d) UX Ori-type circumstellar extinction 
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Fig. 6. M-band light curve, u — r light curve, and rvs. u — r variations for four different objects in NGC 2264: a) CSIMon-000877, a representative of 
the stochastic accretion bursters class; b) CSIMon-001247, a prototypical WTTS; c) CSIMon-000660, an AA Tau-like CTTS; d) CSIMon-000433, 
a UX Ori-like object. Light curves report all flux measurements obtained from the monitoring campaign for a better overview of the variability 
pattern, but only observing sequences obtained in photometric conditions (see Table[Tll are retained for the color analysis. Each point in the r vs. 
u — r diagram for a given object corresponds to the average u — r color and r-band magnitude measured from a single observing sequence; error 
bars are associated with each point based on the lowest rms level detected for field stars at the corresponding brightness (see Fig.[2l. The effects of 
interstellar reddening are traced as a vector in the bottom left quarter of each (r, u — r) diagram. 


dening vector, which is the same in each plot). A more com¬ 
plex trend i s observed for case c), which illustrates an AA Tau- 
like CTTS (lAlencar et al.ll201^ iMcGinnis et al]|2015h . that is, 
an object whose main variability features are driven by quasi- 
periodic occultations of part of the stellar surface from an inner 
disk warp. The r vs. m — r diagram reveals two separate contri¬ 
butions of different nature to the overall variability profile: one 
with color signatures as expected for accretion spots (segment 
at linear increase of u — r with r, with a slope consistent with 


that shown in panel a)), and the other displaying gray brightness 
variations (segment flat in m — r), which is indicative of opaque 
material that crosses the line of sight to the source. These behav¬ 
ior changes occur within only two weeks and average to a global 
ArlA{u — r) slope consistent with that traced by the reddening 
law. Another interesting behavior is shown in panel d). As in 
the previous case, the variability of this system is dominated by 
circumstellar extinction. The system generally becomes redder 
when fainter, with the exception of the first segment of the light 
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Fig. 7. Comparison of photometric amplitudes and Ar/A{u — r) slopes 
observed for representative objects of three different types of NGC 2264 
members (stochastic accretors, WTTS, and AATau-like stars) whose 
variability is driven by accretion bursts (red lines), cold magnetic spots 
(blue lines), or disk occultation (green lines), respectively. Luminosity 
increases from the bottom up along the y-axis and objects become bluer 
from right to left along the x-axis. 

curve, where a decrease in luminosity down to the minimum de¬ 
tected is accompanied by a marked decrease in m — r (i.e., bluer 
colors). Again, two different components can be distinguished 
on the r vs. u — r diagram. At first, the system becomes redder as 
it proceeds from maximum to a mean brightness state, with an 
overall trend consistent with what is expected for an accreting 
source; however, as the brightness drops below a certain level, 
the system is suddenly found to depart from the previous color 
trend and bec ome bluer. This effe ct is reminiscent of the UX Ori 
phenomenon (iHerbst et al.lll994) . 

Individual amplitudes and slopes may vary broadly among 
cases that share the same physical origin (as illustrated in Ap¬ 
pendix]^ for a subsample of WTTS). Nevertheless, we would 
expect to observe noticeable differences between typical color- 
magnitude trends for distinct variable grou ps. Thus, we selecte d 
good representatives of the WTTS, burster (iStauffer et al.ll2014l) . 
and AA Tau-like (iMcGinnis et al.ll2015h classes. This selection 
favored objects in each group that showed the largest photomet¬ 
ric amplitudes in CFHT bands, or equivalently were the least 
affected by noise and hence provided the clearest depiction of 
characteristic variability features. Amplitudes in color and in 
magnitude for these objects are compared in Fig.]?] WTTS show 
the smallest amplitudes and remarkably little color variations as 
the system moves from the brightest to the faintest state. Bursters 
show magnitude variations that are similar to those displayed 
by WTTS in the r band, but are associated with significantly 
larger color variations. The strongest variability amplitudes are 
observed for AA Tau objects, and the average Ar/A(u — r) slope 
detected on a timescale of weeks is found to be consistent with 
the slope expected for interstellar extinction. 

These relative trends are confirmed when analyzing the color 
slopes for the whole sample of variable NGC 2264 members, 




Ar/A(u-r) 


Fig. 8. Histograms of Ar/A{u — r) amplitude ratios are compared for 
WTTS (blue), all CTTS (magenta), accretion-burst-dominated CTTS 
(orange), and circumstellar extinction-dominated CTTS (green). His¬ 
togram bars for bursters and the variable extinction group are shifted by 
—0.025 and 0.025, respectively, relative to the bin center and are placed 
side by side for a better visual comparison of the overall distributions. 
The slope predicted by the interstellar extinction law is shown as a black 
dashed line. 


Table 3. Average color slopes for different typologies of light curves. 


Object class 

No.-' 

Median . 

IQR2 

WTTS 

207 

1.0 

0.7 

CTTS 

212 

0.6 

0.5 

U. Bursters 

18 

0.5 

0.4 

U. Circumstellar extinction 

26 

0.9 

0.5 


Notes. For comparison, the inter stellar extinction slope is 1.22 (from 
data in iFiorucci & Munarill2003h . Number of objects in the corre¬ 
sponding class. Interquartile range. 


WTTS or CTTS, with a particular focus among the latter on ob¬ 
jects showing either a stochastic burster or an AA Tau-like na¬ 
ture. Median slope and statistical dispersion measured in each of 
these four groups are listed in Table |3] (while individual values 
of Ar/A{u — r) are reported in Table |2] and corresponding dis¬ 
tributions are illustrated in Fig.j^. CTTS show on average shal¬ 
lower slopes than WTTS; this difference is more marked when 
only considering accretion-dominated objects. AA Tau objects 
and WTTS display similar slope ranges, in turn consistent with 
the interstellar extinction slope. 

3.4. Spot models 

As di scussed by several authors over the past few year s (e.g.. Pa¬ 
per F iCostigan et al.ll2()T4l iGrankin et al.ll2007il200^ . baselines 
of days to weeks appear to be the leading timescales for YSO 
variability (at least up to several years). A major contribution 
to this variability arises, on average, from rotational modulation 
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due to surface inhomogeneities, whose nature is linked with the 
magnetic activity of the star and/or with ongoing disk accretion. 

Interesting indications of the properties and dominating fea¬ 
tures of individual objects can thus be inferred by attempting to 
reproduce the observed multiwavelength variability signatures 
with spot models. In this picture, the variability is assumed to 
arise from a region (single spot or spot distribution) of different 
temperature at the stellar surface, whose emission modulates the 
observed luminosity of the star while it rotates. 

Several studies, including extensive photometric monitoring 
and spectropolarim etric observations of individual objects (e.g., 
iDonati et al.ll20ir)h . have shown that the surface of TTS is pos¬ 
sibly covered by multiple spots or spot groups. While our ca¬ 
pability of inferring a detailed description of surface spot distri¬ 
bution is limited by the poor geometric constraints, useful indi¬ 
cations on effective temperature and area covered by spots can 
be deduced from the observed variability amplitudes at differ¬ 
ent wavelengths. As mentioned earlier, both cold (magnetic ac¬ 
tivity) and hot (accretion) spots at the stellar surface will pro¬ 
duce the same qualitative effects; objects are expected to be red¬ 
der at fainter states and variability amplitudes will be larger at 
shorter wavelengths. Spot parameters (notably temperature and 
fractional area coverage relative to stellar photosphere) are, then, 
uniquely determined by differential color variations, that is, by 
the rate at wh ich the observed a mplitudes decrease toward longer 
wavelengths (IVrba et al.lll993h . 


3.4.1. Model adopted in this study 


To explore the nature of variability across the PMS population 
of NGC 2264, we followed the approach of iBouvier et al.l (Il993h 
and adopted a spot model that does not introduce any assump¬ 
tions of the number and shape of spots at the stellar surface. The 
model assumes, however, i) a uniform temperature for all spots 
at the stellar surface, and ii) a blackbody distribution for stellar 
and spot emission. 

As detailed in Bouyier et al.l (Il993h . after IVogd (Il98lh and 
[Torres & Ferraz Mellol ( 1973h . modulated variability amplitudes 
due to surface spots can be described as 


Am{A) 


-2.5 log 



Teq 


’ 


(3) 


where S(/l) is the specific intensity of the immaculate photo¬ 
sphere at the stellar disk center, S/(/l) is the same quantity for 
the spotted region, 5'(d)/5(/l) = B;i{Tspot)/B,x{Tphot) in the 
blackbody approximation, Teq is a lower limit to the true frac¬ 
tion of stellar surface covered by spots, and ^{A) is the linear 
limb-darkening coefficient. 

Photospheric temperatures Tphot were inferre d from the spec¬ 
tral ty pes of the objects following the scale of ICohen & Kuhil 
(Il979ln . while linear limb-darkening coefficients were deduced 


‘ The scale of ICohen & Kuhil dl979li is similar to that of 

iKenvo 

n & Hartmann 

Jl99^ for spectral types earlier than Ml and to 
(l2003b for later spectral types. It differs somewhat 

that of 

Luhman et al. 


from the more recent determination of iPecaut & MamaiekI l l2013h for 
PMS stars, notably in th e earlier (i.e., K ) spectral ty pe range. Our choice 
of favoring the scale of ICohen & Kuhil over that of iPecaut & MamaiekI 
here is motivated by the fact that the color sequences on which the for¬ 
mer is based provide a better match to the empirical cluster sequence in 
CFHT color-magnitude diagrams, especially for earlier type stars, than 
models relevant to the latter. 


from the compilation of IClaret & Bloemi^ (l201 ih . assuming so¬ 
lar metallicity and log g-A. Limb-darkening coefficients tabu¬ 
lated for spectral type M2 were uniformly extended to later spec¬ 
tral types that were not sampled in the reference compilation at 
the metallicity and gravity values adopted. Spot parameters Tspot 
and Teq are thus the only unknowns in Eq.|3]for a given object 
and were deduced from the model by fitting the expression for 
Am to the observed u- and r-band variability amplitudes. 

The limitations of this model and t he implication s of the as¬ 
sumptions are discussed extensively in iBouvier et al.l (Il993l) . 

To find the (Tspot, f^eq) pair that best describes the variabil¬ 
ity properties observed for a given object, we explored a two- 
dimensional grid, ranging from lOOOK to lOOOOK in T^pot and 
from 0.1% to 90% in Teq, with a step of 5 K in Tspot and of 
0.002% in 'Feq- This parameter space range was enlarged when 
the best solution fell at the edge of the explored domain. For 
each point of the grid, the theoretical amplitudes Am defined as 
in Eq.[3]were computed, and their agreement with the observed 
variability amplitudes was estimated through statistics; the 
best solution is defined as the one that minimizes 

The significance of the result was evaluated by investigat¬ 
ing how stable this “best” solution is relative to weak variations 
of the input “observed” amplitudes. Spot model predictions are 
typically affected by a de gree of degeneryy bet ween different 
spot parameters (see, e.g.. lWalkowicz et'^l2013l) : a given set of 
amplitudes could in principle be reproduced either in terms of 
more extended surface spots with a modest temperature differ¬ 
ence relative to the stellar photosphere, or in terms of a smaller 
spot distribution with a larger temperature difference. To esti¬ 
mate the uncertainty on the spot model result for a given object, 
we did the following: 

1. we took the u- and r-band amplitudes (amp,„ amp^) and the 
associated error bars err„, err^. (defined as err,„ = y/2o^, 
where cr„, is the photometric noise at the relevant brightness 
location, traced as a fit to the lower envelope of the appropri¬ 
ate field stars distribution in Fig.|2]i; 

2. we used the rejection method for generating rando m devi¬ 
ates within a normal distribution dPress et al.lll992h to pro¬ 
duce 200 random (amp[f^', amp'^'^’O pairs, where amp^® ” is ex¬ 
tracted from a Gaussian distribution centered on amp„, with 
standard deviation err,„; 

3. for each of the randomly generated (m, r) amplitude pairs, we 
ran the spot model in Eq.[3and computed the spot parameters 
that best reproduced these test variability amplitudes; 

4. the spot model results for each of the random amplitudes 
pairs were used to build a Tspot and a Teq distribution, whose 
mean and standard deviation were extracted as a best value 
and a corresponding uncertainty, respectively. 

3.4.2. Global picture of TTS spot properties in NGC 2264 

An overview of the different nature of modulated variability ex¬ 
hibited by accreting and no n-accreting members of NGC 2264 is 
provided in Fig.|9](see also lBouvier et alJll995h . This compares 
the u- and r-band amplitudes for different classes of members to 
cool (upper panel) and hot (lower panel) spot model predictions, 
computed assuming, for illustration purposes, a typical photo- 
spheric temperature of 4000 K and varying spot temperatures 
and filling factors. 

A first interesting feature of these two diagrams is that u- 
band variability amplitudes are always larger than r-band ampli¬ 
tudes. This is consistent with expectations in the spot description 
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2.5 - Cool-Spot models 
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2.5 - Hot-Spot models 
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Fig. 9. u- and r-band variability amplitudes observed for WTTS (blue squares) and CTTS (magenta circles) are compared to cool (upper panel) 
and hot (lower panel) spot model predictions. Model tracks are traced assuming a fixed photospheric temperature Tpi,„, of 4000 K (corresponding 
to a K7 star). Each model track corresponds to a different spot temperature Tjpn,: in the upper panel (cool spots), four different Tsp„t values are 
considered, ranging from 200 K to 1500 K lower than F,in the lower panel (hot spots), five different T^po, values are explored, ranging from 
200 K to 4000 K above the photospheric temperature. Along each track, varies from 0.1% to 60%; reference values of the spot-to-star surface 
ratio are marked wit h crosses for clarity. For CTTS, two additional subclasses are h ighlighted: objects do minated by stochastic accretion bursts 
dStauffer et al.ll2014l orange dots) and objects dominated by circumstellar extinction dMcCinnis et alj2015l . green triangles), either quasi-periodic 
(AA Tau objects) or aperiodic. 


of modulated variability. WTTS typically display lower variabil¬ 
ity amplitudes than CTTS, < 0.4 mag in r and 0.5 mag in u. Cool 
spot model tracks, depicted in the upper panel of Fig.0 span the 
WTTS locus over the entire amplitude range fairly well; this sug¬ 
gests that in most cases, WTTS variability can be convincingly 
reproduced in terms of cool magnetic spots, having a tempera¬ 
ture of ~500-1000 degrees lower than Tpi,ot and covering up to 
~30% of the stellar surface. A smaller fraction of WTTS are 


located in the region of the diagram dominated by hot spots; 
enhanced chromospheric activity may explain the photometric 
properties observed for some of these objects, while others might 
actually be accreting stars at levels too low to be detected. A 
non-negligible overlap exists between the distributions in am¬ 
plitudes of WTTS and CTTS; however, a more significant frac¬ 
tion of CTTS is observed to exhibit larger variability amplitudes 
than WTTS, and, at a given value of amp, , the bulk of accreting 
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^phot ^spot 

Fig. 10. Best spot model solutions for CTTS (magenta circles) and 
WTTS (blue squares) in NGC 2264 are statistically compared. A typ¬ 
ical error bar is shown as a black cross in the upper left comer of the 
diagram. Error estimates result from the uncertainties on the photomet¬ 
ric amplitudes, following the procedure described in Sect. l3.4.T] plus a 
lower order correction accounting for uncertainties on the photospheric 
temperature (Tpi,o,) and on the shape of the limb-darkening coefficients 
^(/i). The red vertical line marks the zero-value temperature difference 
between stellar photosphere and spot distribution. For CTTS, objects 
dominated by stochastic accretion bursts are marked in orange, while 
green triangles correspond to AA Tau-like objects. Histograms shown 
in side panels at the top and to the right of the diagram compare the 
CTTS (magenta) vs. WTTS (blue) distribution in spot properties (tem¬ 
perature difference relative to the stellar photosphere and percentage of 
surface covered, respectively). The heights of histogram bars are nor¬ 
malized to the highest channel count in each distribution. 


members is located at larger M-band amplitudes than their non¬ 
accreting counterparts. Cool spot model predictions are not able 
to reproduce the color variability of objects located in this part 
of the CTTS distribution; a better fit to their photometric am¬ 
plitudes is provided by hot-spot model tracks in the lower panel 
of Fig.|9] Temperature differences between hot spots and stellar 
photosphere range from <1000 K to several thousand K. Typical 
hot-spot distributions extend over ~5-10% of the stellar surface 
for spots 2000 K-hotter than the photosphere and down to a few 
percent for the hottest spots. 


3.4.3. Spot model results 

Figure[T0] synthesizes the main features of the spot properties 
inferred, case by case, for the populations of CTTS and WTTS 
in NGC 2264, using the model described in Sect. l3.4.TI 

Three distinct loci can be identified in this diagram. The first, 
located to the right of the red line (i.e., at positive values of 
Tphot- Tspot), below =40%, is the main locus populated by 
WTTS. In most cases, the best spot solution suggests a surface 
spot distribution ~500 K colder than the stellar photosphere, ex¬ 
tending over ~ 10-30% of the stellar surface. A small number of 
CTTS are also located in this region of the diagram, with prop¬ 


erties consistent with those deduced for WTTS, which suggests 
that cool spots of magnetic origin are still a primary cause of 
modulated variability in these cases. 

The second locus is located to the left of the red line (i.e., at 
negative values of Tphot - Tspot) and is the locus primarily pop¬ 
ulated by CTTS. Spot distributions that best reproduce the ob¬ 
served photometric amplitudes in individual cases are typically 
~ 1000- 2000 K hotter than the stellar photosphere, with a tail of 
objects for which temperature differences of over 3000-4000 K 
are predicted by spot model results. Although the Tspot - Teq de¬ 
generacy in spot model results may play a non-negligible role in 
the anticorrelation between these two quantities, observed both 
to the left and right of the red line in Fig.[T0] a significant differ¬ 
ence between the distributions of CTTS and WTTS over Teq can 
be detected, with the former being largely located below 5-10%, 
while the latter are mainly comprised between 5-10% and 30%, 
as discussed earlier. 

The third locus in the diagram is the point distribution around 
the Tphot - Tspot - 0 line, above - 40%. This is mainly popu¬ 
lated by CTTS. These objects stand out for the small tempera¬ 
ture difference between photosphere and best spot model and for 
the corresponding high value of significantly higher than 
the typical values for WTTS or CTTS. This suggests that for this 
subgroup of objects, spot models are ill-suited to describe the ob¬ 
served amplitudes and color properties, and hence other compo¬ 
nents may be dominating the observed variability. This conclu¬ 
sion is supported by several AA Tau-like objects in this group, 
whose variability is known to be dominated by quasi-periodic 
disk occultation. 


4. Discussion 

4.1. Photometric behavior of CTTS and WTTS at short 
wavelengths 

From a statistical point of view, accreting YSOs (CTTS) exhibit 
significantly different variability properties from disk-free young 
stars (WTTS). Higher levels of variability are associated with 
the former group, both in the optical (r band) and, even more 
markedly, at UV wavelengths (u band); typical photometric am¬ 
plitudes measured across the CTTS population amount to about 
three times those characteristic of the WTTS group. 

A definite anticorrelation trend is observed between variabil¬ 
ity amplitudes and stellar mass for WTTS, both in the r and in the 
u band, as illustrated in Fig.[TT] This would suggest that smaller 
and/or more uniformly distributed spots, translating into more 
modest peak-to-peak luminosity variations, are found at the sur¬ 
face of more massive PMS stars than for lower mass objects. 
Conversely, no significant correlation properties are observed be¬ 
tween mass and photometric amplitudes measured for CTTS. 

Bluer CTTS (i.e., those displaying higher UV excess rela¬ 
tive to photospheric colors) statistically tend to exhibit stronger 
UV variability. The UV excess is one of the most direct diag¬ 
nostics of the accretion process; this statistical correlation, and 
the different variability behavior pertaining to the CTTS popu¬ 
lation compared to WTTS, suggest that the variability scenario 
observed for CTTS is dominated by an active circumstellar disk 
and not by magnetic activity common to non-accreting young 
stars (WTTS). 

Variability amplitudes measured for WTTS are best inter¬ 
preted in terms of geometric luminosity modulation from cool 
surface spots, a result of the magnetic activity of the star. Typi¬ 
cal spot parameters found in this study for WTTS amount to an 
effective temperature ~500 degrees colder than the photosphere, 
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Fig. 11. Photometric amplitudes in the u band (upper panel) and in the r 
band (lower panel) as a function of stellar mass are compared for CTTS 
(magenta circles) and WTTS (blue squares). Mass estimates are derived 
as described in Paper I. 


and area coverage from 5-10% to 30% of the stellar surface. 
Conversely, CTTS variability appears to be dominated by hot 
accretion spots, with typical effective temperatures up to 4000 K 
hotter than the stellar photosphere, and that extend over a few 
percent of the stellar surface. 

It is important to point out that although CTTS and WTTS 
appear to be statistically distinct in many metrics, the two pop¬ 
ulations display a non-negligible overlap in photometric behav¬ 
ior. Some CTTS, classified as such based on “standard” indica¬ 
tors of accretion disk (H„ emission, near-IR or UV excess), are 
found to exhibit photometric properties consistent with objects 
that lack disk signatures, and vice versa. This would indicate that 
a continuum of possible behaviors instead of a clear-cut transi¬ 
tion exists between the state of a “classical” T Tauri and that of 
a diskless, weak-lined T Tauri star. Figure 16 of Paper I clearly 
shows that a broad variety of accretion regimes coexists within 
the same young stellar population (from strong accretion rates 

Mace ~ lO^^MQ/yr down to marginal Mflcc ^ 10~*°MQ/yr). On 

the other hand, some sources classified as WTTS by our criteria 
might actually be objects with low undetected accretion or in a 
“quiescent” state. Cases of WTTS-like young stars that episodi¬ 
cally cross the C TTS-defining borders have indeed been reported 
in the literature (iMurphv et al.ll20llllCieza et al.ll2013h . 

4.2. Color trends 

Independently of their specific nature, monitored members are 
typically found to become redder as they fade, and observed 
variability amplitudes are larger at shorter wavelengths. This 
property is well understo od when considering a spot-dorninated 
variability sc enario (e.g., iBouvier et al.lll993[ IVrba et aDll993l: 
iHerbst et al.lll994 . 

While this qualitative trend is common to systems of dif¬ 
ferent nature, the rate at which the amplitudes decrease with A 


is more informative regarding the specific mechanisms at play. 
Different behaviors are indeed observed when comparing the 
typical Ar/A(M — r) slopes measured for the WTTS vs. CTTS 
classes. Accretion-dominated objects stand out for their signifi¬ 
cantly shallower slopes, which are indicative of a much stronger 
contrast at UV wavelengths than that measured at longer wave¬ 
lengths. Ambiguous results are obtained for the WTTS class and 
the subgroup of disk-bearing objects (AA Tau objects), whose 
variability appears to be dominated by circumstellar extinction. 
Typical slopes measured in both cases are found to be consistent 
with each other and with the color slope predicted for interstel¬ 
lar extinction; a somewhat larger dispersion in values is detected 
across WTTS. A similarity in color behavior between what is ob¬ 
served for WTTS and what is expected in the case of circumstel¬ 
lar extinction (assum ing interstellar extinc tion properties) was 
previously noted by iGrankin et al.l (l2008h in the optical. We 
stress that while a linear trend appears to provide a detailed de¬ 
scription of the color variation with magnit ude for spotted stars 
(as al so observed, on longer timescales, bv IGrankin et alJl2007i 
l2008l) . this merely provides average information on the actual 
color behavior observed for some AA Tau-like or circumstel¬ 
lar extinction-dominated systems (as illustrated in panels c) and 
d) of Fig. 13. Indeed, complex and ph ase-varying color beh avior 
is well documented in AA Tau itself (iBouvier et al.ll20()3h . The 
alternation of phases of colored and gray magnitude variations 
for these objects may indicate non-uniform extinction properties 
across the occulting screen (e.g., more opaque at the center and 
optically thin at the edges). Furthermore, increased veiling at the 
epochs of maximum accretion shock visibility (i.e., close to the 
occultation event, in the assumption that the inner disk warp is 
located at the base of the accretion column) may determine a 
sudden transitio n of the system betwe en different color regimes. 

Remarkably, IGrankin et al.l (l2007h detected color slopes for 
CTTS that are quit e similar to those measured for WTTS 
(IGrankin et al.l1200^ in optical bands {V, R). This result, com¬ 
pared with the present analysis, shows that hot and cold spots are 
primarily distinguished at short wavelengths. 

4.3. Stability and leading timescales foru-band variability 

The present study focuses primarily on the exploration of UV 
variability on timescales relevant to stellar rotation. However, 
comparing variability signatures on different timescales is instru¬ 
mental in achieving a more complete picture of YSO dynamics. 

The short timescales (hours) show little variability com¬ 
pared to the amount measured on days/weeks (more mid-term 
timescales). To estimate a typical amount of short-term vari¬ 
ability, we computed the average rms scatter of photometric 
measurements obtained at different epochs within single nights 
(Sect.|2]l across our sample. CTTS are typically found to dis¬ 
play an rms variability of 0.06 mag on these time baselines, cor¬ 
responding to barely 8-9% of the amplitudes measured during 
the full CFHT monitoring. Similar results are obtained for the 
WTTS group, with typical short-term rms variability of 0.025 
mag, or about 10% with respect to the corresponding mid-term 
variability amplitudes. 

This result was partly introduced in Sect. 4.3 of Paper I in the 
context of a more quantitative assessment of the contribution of 
modulated variability to the total amount of variability detected 
for accreting objects. We found that up to ~75% of the UV 
variability measured for typical CTTS on week-long timescales 
(generally a few to several tenths of magnitude) is simply due 
to geometric effects of rotational modulation. Of the remaining 
fraction linked with intrinsic variability, a larger contribution can 
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Fig. 12. n-magnitude differences between single-epoch photometry in 
December 2010 and median light curve photometry in February 2012 
are compared to the semi-amplitudes of u-band light curves from Febru¬ 
ary 2012 monitoring for WTTS (blue squares), CTTS (magenta circles), 
and CTTS accretion-dominated (orange dots) and extinction-dominated 
(green triangles) subgroups. The equality line is traced in black to guide 
the eye. 

be attributed to cycle-to-cycle variations (possibly related to ac¬ 
cretion spot evolution), while a smaller component is statistically 
contributed on hour-long timescales, which is more sensitive to 
unstable, stochastic behaviors in the accretion process. 

To evaluate the effect of mid-term timescales on long¬ 
term (years) variations, we compared u-band photometry ob¬ 
tained during the CFHT/MegaCam monitoring survey of Febru¬ 
ary 2012 with single-epoch photometry obtained with the same 
instrument during a preliminary, single-epoch mapping survey 
performed in December 2010. Details on this first photometric 
run and on subsequent data reduction are provided in Paper I. 
If mid-term timescales dominate the variability pattern exhibited 
by TTS on longer baselines, we expect the photometric mea¬ 
sure from the snapshot survey of December 2010 (m 2010 ) to cor¬ 
respond to a random epoch in the light curve of the object re¬ 
constructed from 2012 monitoring. Hence, when measuring the 
absolute difference between u 2010 and the median light-curve 
magnitude from February 2012 expect this differ¬ 

ence to be smaller than (or similar to) half the variability am¬ 
plitude amp u of the light curve. Conversely, if other variability 
mechanisms are predominant on year-long timescales, we ex¬ 
pect to measure 1 1/2010 — “ 2012 ! values on average higher than 
amp u/2. This comparison is shown in Fig. [12] for all variable 
members considered in this study. A clear offset along the equal¬ 
ity line is observed between the WTTS and the CTTS groups in 
the diagram, with the latter group located at signihcantly higher 
variability amplitude values. However, there are no signihcant 
differences between the two stellar groups in the indicator of 
interest here, that is, the |m 2 oio ^ ratio. Both 

groups are distributed in clouds of similar properties about the 
equality line on the diagram, with a mean ratio of 0.8 (i.e., typi¬ 
cal epoch-to-epoch difference that can be fully accounted for by 


mid-term variability) and comparable rms scatter. This suggests 
that similar timescales or mechanisms (notably rotational mod¬ 
ulation) are responsible for the long-term variations observed in 
the two cases. 

For the CTTS group, we additionally probed the stability of 
M-band variability on longer timescales (>10 years) by compar¬ 
ing the UV excess luminosity measured and monitored in our 
survey (see Paper I) with values measured in a si milar single¬ 
epoch U -band survey of NGC 2264 performed by iRebull et al.l 
(l2002h . Results of this comparison are shown in Fig. 13 of Pa¬ 
per I. Again, we observe that median UV exces s luminositie s 
from 2012 and single-epoch measurements from iRebull et all ’s 
study distribute around the equality line with an rms scatter 
consistent with the typical amount of variability detected on a 
timescale of weeks. 

These results suggest that the mid-term timescale is the lead¬ 
ing timescale for YSO variability, at least up to baselines of sev¬ 
eral years. The similarity of WTTS and CTTS in this respect also 
suggests that the typical dynamics of the accretion process on 
CTTS is intrinsically stable over timescales of years, although 
it may be variable in the shorter term. In other words, single 
accretion events may be shorter lived, and the dynamics of in¬ 
teraction between disk and stellar magnetosphere may lead to 
rapidly evolving surface spots (on timescales of a few rotational 
cycles or less) and irregular light curve morphology. However, 
the averaged spot properties, indicative of the intensity of mass 
accretion onto the star and of the mechanisms that regulate the 
process, can persis t over timescales of y ears. Similar conclusions 
were presented bv iGrankin et alJ (l2007h . 


5. Conclusions 

The work reported here provides a hrst extensive mapping of 
variability properties at short (UV) wavelengths for a whole star¬ 
forming region (the NGC 2264 cluster) and its several hundred 
members. The sample encompasses over 750 young stars in the 
region, of which about 40% show evidence of disk accretion. 
Variability was monitored simultaneously in the optical (r band) 
and in the UV (« band) over two full weeks, with multiple mea¬ 
surements per observing night. Our statistical analysis globally 
confirms and strengthens earlier results on the manifold nature of 
variability displayed by young stellar objects. Disk-bearing clus¬ 
ter members (CTTS) exhibit significantly higher levels of vari¬ 
ability than non-accreting objects (WTTS) on timescales rang¬ 
ing from hours to weeks. A statistical correlation is observed 
between the amount of photometric variability and that of UV 
excess measured across the former group, which indicates that 
ongoing accretion and star-disk interaction are the main driv¬ 
ing factors of the variable nature of CTTS. Amplitude ranges 
and color variations monitored for WTTS are best reproduced in 
terms of cool surface spots linked with magnetic activity, which 
modulate the luminosity of the star as it rotates. Hot-spot mod¬ 
els are instead required to account for the stronger variability and 
steeper amplitude increase from the optical to the UV observed 
for typical CTTS. In both groups, similar amounts of variabil¬ 
ity are observed in the mid- (weeks) and long term (years); this 
indicates that rotational modulation is the main source of YSO 
variability over several hundred periods, independently of their 
accretion status. This in turn suggests that the underlying phys¬ 
ical processes (namely magnetic activity in WTTS and disk ac¬ 
cretion in CTTS) have a relatively stable nature in the longer 
term, although disk accretion can exhibit more erratic behavior 
in the shorter term. 
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Appendix A: Spotted stars and color slopes 

Individual amplitudes and color slopes may vary broadly among 
cases that share the same physical origin for the observed vari¬ 
ability; this reflects intrinsic object-to-object differences in the 
relative surface extent and temperature contrast provided by the 
intervening layer. This point is illustrated in Fig. IA.il for a sub¬ 
set of WTTS members from our sample with uniform spec¬ 
tral type of Ml (hence homogeneous in photospheric tempera¬ 
ture). For these diskless young stars, the observed variability is 
most straightforwardly interpreted in terms of surface magnetic 
spots, cooler than the stellar photosphere, that modulate the pho¬ 
tospheric emission. Figure lATI shows that larger r-band ampli¬ 
tudes of variability typically correspond to steeper Ar / A(m — r) 
slopes; this suggests that a less significant T-dependence in the 
spot-to-star contrast is observed for objects with the largest vari¬ 
ability amplitudes. This result can be understood if we con¬ 
sider that along the average trend, larger amplitudes reflect larger 
temperature differences between spot and stellar photosphere. 
Darker spots emit less flux, resulting in more marked luminos¬ 
ity differences between the phases of minimum and maximum 
spot visibility during stellar rotation. This luminosity difference, 
however, will be less colored, because the emission spectrum for 
darker spots peaks at longer wavelengths, and hence progres¬ 
sively contributes less in the observed d-window; consequently, 
the net effect of darker spots in the spectral range of interest will 
be a more uniform flux decrease when the spot is in view. Simi¬ 
lar trends to that shown in Fig. lA. H are found at different spectral 
types; the scatter in amplitudes around the average trend may 
reflect a variety in the effective percentage of stellar surface cov¬ 
ered by spots (for a given temperature, spot distributions cover¬ 
ing a larger area will produce larger variability amplitudes, but 
the color dependence of the contrast will remain unchanged). In 
addition, some marginal accretion activity might still be present 
in the more variable stars here classified as WTTS, hence con¬ 
tributing to the vertical scatter in values. 



Fig. A.l. Measured r-band amplitudes (Ar) and corresponding slopes 
(Ar/A{u — r)) are compared for WTTS members of spectral type Ml in 
our sample. Error bars on photometric amplitudes follow from the pho¬ 
tometric noise estimate based on field stars sequence in Fig.|2l error bars 
on slope values are computed by simulating 1 000 random (Ar, A(m —r)) 
pairs uniformly distributed around the observed values within the uncer¬ 
tainties and measuring the corresponding scatter on computed slopes. 
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